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ABSTRACT 



Context. Planet perturbations have been often invoked as a potential explanation for many spatial structures that have been imaged 
in debris discs. So far this issue has been mostly investigated with pure N-body numerical models, which neglect the crucial effect 
collisions within the disc can have on the disc's response to dynamical perturbations. 

Aims. We numerically investigate how the coupled effect of collisions and radiation pressure can affect the formation and survival of 
radial and azimutal structures in a disc perturbed by a planet. We consider two different set-ups: a planet embedded within an extended 
disc and a planet exterior to an inner debris ring. One important issue we want to address is under which conditions a planet's signature 
can be observable in a collisionally active disc. 

Methods. We use the DyCoSS code of Thebault(2012), which is designed to investigate the structure of perturbed debris discs at 
dynamical and collisional steady-state, and derive synthetic images of the system in scattered light. The planet's mass and orbit, as 
well as the disc's collisional activity (parameterized by its average vertical optical depth t ) are explored as free parameters. 
Results. We find that collisions always significantly damp planet-induced spatial structures. For the case of an embedded planet, the 
planet's signature, mostly a density gap around its radial position, should remain detectable in head-on images if M plalu , t > M Satum . If 
the system is seen edge-on, however, inferring the presence of the planet is much more difficult, as only weak asymmetries remain in 
a collisionally active disc, although some planet-induced signatures might be observable under very favourable conditions. 
For the case of an inner ring and an external planet, planetary perturbations cannot prevent collision-produced small fragments from 
populating the regions beyond the ring. The radial luminosity profile exterior to the ring is in most cases close to the one it should 
have in the absence of the external planet. The most significant signature left by a Jovian planet on a circular orbit are precessing 
azimutal structures that can be used to indirectly infer its presence. For a planet on an eccentric orbit, we show that the ring becomes 
elliptic and that the well known pericentre glow effect is visible despite of collisions and radiation pressure, but that detecting such 
features in observed discs is not an unambiguous indicator of the presence of an outer planet. 

Key words, stars: circumstellar matter - planetary systems: formation - stars: individiual (HR4796, Fomalhaut, HD202628) 



1. Introduction 

Most imaged dusty debris disc£] display pronounced radial 
and azimutal structures, which are the sign that these sys- 
tems are dynamically active and that "something" is shaping 
them. Depending on the type of spatial features, i.e ., two-sided 
asymmetries, spirals, warps, c l umps or rings (e.g [Kalas et al.l 
I2005L iGolimovski elaH l2006t ISchneider et all I2009-I) . several 
possibl e explanations have been i nvestigated: transient v iolent 
events (Kenyon & Bromley 2005; Grigori eya et al.l l2007). cou- 
pling to gas drag dTakeuchi & Artvmowiczll200ll). companion 



star perturbations ( Augereau & Papaloizou 2004; Thebault et al 
2010; Th ebaultl2012l) or dynamically cold discs dThebault & Wu 



depending on resonance order and pl anet orbits, pronounced az- 
imutal and radial overdensities (e.g. iKuchner & Holmanl [2003; 
IWvatll 120061; iReche et al] 120081) . Likewise, the chaotic region 
surrounding an embedded planet can efficiently truncate a disc, 
inducing sharp inner or outer edges and ring-like structures 
dWisdoml fl980; M ustill & Wvattll2012l) . Embedded planets can 
also trigger transient spiral structures that can often be long- 
lived enough to be observed in extended debris discs dWvattl 
2005). Finally, inclined planets have been found to be a possi- 
ble cause for warped discs s uch as in /?-Pictoris ( Mouille t et al 
1997]: lAugereau etaT]l200lb iDawson et aDl201 lb IChauvin et al 



2012) 



2008). However, the most commonly proposed scenario for the 



vast majority of theses structures is the presence of (usually 
unseen) perturbing planets. The dynamical effect of a(several) 
planet(s) on a disc has been investigated in numerous studies, 
usually numerical investigations based on deterministic N-body 
codes. These investigations have lead to several important re- 
sults, notably that dust can easily be captured in mean motion 
resonances with inner planets that are migrating, thus creating, 
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However, caution is required when trying to directly match 
N-body code results to imaged discs. Indeed, most resolved im- 
ages are obtained in scattered light, at wavelengths for which 
the flux is dominated by the smallest dust grains on bound or- 
bits, which are strongly affe cted by stellar radiation pressure 
dThebault & Augereau! 12007). While including the effect of ra- 
diation pressure in N-body codes is not difficult in itself, usually 
by correcting the stellar gravity by a factor (1 - ft) 05 (where 
/3 is the ratio between radiation pressure and stellar gravity), 
the fact that it is a size-dependent effect requires to have an 
idea of the size distribution in the disc. The problem is that 
the size distribution is imposed by the collisional evolution of 
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the system, and that collisions are not taken into account in 
pure N-body codes. To circumvent this problem, an idea would 
be to perform series of dynamics+radiation-pressure N-body 
runs, each for a given particle size, and then recombine them 
with a weight given by theoretical collisional-equilibrium size 
distributions; for exam ple th e classical differential distribution 
in -3.5 of rPohnanvil <ll969h (iMoro-Martm & Malhotral [ 2002: 
lOuillen & Thorndikdl2002l; I Wvattl 120061: lErtel et al.ll2012h . But 
the implicit assumption behind this procedure, i.e., that size 
distributions are not affected by spatial structures or dynami- 
cal p erturbations, has been pr oven to be faulty. As a mat ter of 
fact. lStrubbe & Chiang! d2006l) and lThebault & Wul (120081) have 
shown that, even in an unperturbed system consisting of a nar- 
row parent body ring, collisional evolution naturally creates a 
disc-integrated overabundance of small grains close to the blow- 
out limit, and thus a strong departure from a "classical" colli- 
sional equilibrium distribution. Size distributions in perturbed 
systems could depart even more from standard power-laws, with 
potentially strong spatial variations Q. It is thus very difficult to 
know in advance the relative contributions of different grain pop- 
ulations to the system's luminosity. 

The non-inclusion of collisions in N-body studies has other 
problematic consequences. An important issue is that of the 
timescales: if the typical collision times are shorter than the dy- 
namical timescales required for spatial structures to form, then 
these structures' development can be affected or even hindered. 
This is why p ure dynamical mo delling, such as the recent in- 
vestigation by Ertel et af] d2012l) unavoidably overestimate the 
level of spatial structures in dusty discs, especially dense ones. 
Another issue is the unavoidable feedback of the collisions on 
the dynamics, as they dissipate energy and produce fragments 
on new orbits. Last but not least, the steady production of small 
collisional fragments, coupled to the effect of radiation pressure 
on them, can inject matter in dynamically unstable regions be- 
fore planets can efficiently remove them. 

Including collisions, especially fragmenting ones as in de- 
bris discs, into a N-body scheme has proven a very arduous 
task. "Brute force" methods, where test particles are effec- 
tively broken into fragm ents that are then dynamically evolved 
(Beauge & Aa rsethl 1 19901) . are in principle the most reliable 
ones, but lead to an exponential increase of particle numbers that 
is very quickly unmanageable. Another alternative is to run pure 
N-body runs to get an idea of the dynamical effect of planets, 
and then use the results, in the form of averaged laws for impact 
probabilities and collisi on velocities, in classica l particle-in-a- 
box collisional codes dKenvon & B romley 2008, and references 
therein). This approach remains however limited to 1-D spa- 
tial resolution and cannot study the formation of fine structures. 
The most promising approach is probably a hybrid model, where 
each particle of the N-body code is a "super particle" (SP) stand- 
ing for a cloud of real bodies sharing a common physical size, 
whose mutual collisions are then treated wit h a particle-in-a- 
box s cheme. The first versions of such a model (Grigoriev a et al.l 
l2007h had a progressive increase of the number of SPs and where 
limited to s hort timescales. How e ver, th e recent "LIDT" code de- 
veloped by ICharnoz & Taillifel (1201 2|) has solved this problem 
but is so far operational only for the simpler case of low velocity 
(i.e., mostly accreting) collisions in proto-planetary discs. 

So far, only two models have managed to incorporate, to 
some extent, collisional effects into N-body codes in the debris- 



2 We note that there is also observational support for wavelength- 
dependent (and therefore grain size-depen dent) structures in debris 
discs fe.g. lSu et al.ll2005l : IWilner et al.ll201 lh 



disc case of high-velocity fragmenting impacts. Both codes are 
comparable in terms of the level of collisional effects they in- 
clude: they do not fully couple collisions and dynamics in a 
self-consistent way but are designed to study systems perturbed 
by one perturber, once a dynamical and collisional steady state 
has been reached. The first one is the Co llisional Grooming 
Algorithm (CGA) of IStark & Kuchneri (120091) . Its principle is to 
first perform "seed runs" of collisionless particles, from which 
streams of successive particle positions are recorded to construct 
density maps. These maps are then used to derive collision de- 
struction probabilities in new seed runs, from which new density 
maps are derived, etc. The process is iterated until convergence 
is reached. This code is limited (so far) to the case of one inter- 
nal perturber on a ci rcular orbifl but has giv en impressive results 
for the Kuiper Belt (Kuc hner & Sta rk 2010). The second code is 
the "DyCoS S" (for "Dynami cs and Collisions at Steady State") 
algorithm of lThebaultl d2012l hereafter TB012). Although being 
significantly less user-friendly in its use and lacking the intrinsic 
coherence of the CGA, it is (in their respective present versions) 
more versatile with respect to the perturber's orbit, which can be 
circular or eccentric, and external to, interior to or embedded in 
the disc. It is also better suited to the study of short timescale 
processes, such as the competing effect of dynamical removal 
and collisional production of small, radiation-pressure-affected 
grains in high-density discs. 

We use here an updated version of DyCoSS, which had ini- 
tially been developped to study circumprimary discs in binaries 
(TB012), adapted to the case of a planetary perturber. The main 
issue we are planning to address is how the collisional activity 
inside a disc affects the dynamical sculpting a planet excerts on 
it, in particular because of the constant flow of small collision 
fragments flying around the system because of radiation pres- 
sure. One crucial point is under which conditions a planet's sig- 
nature can be detectable in a collisionally active disc, and if this 
signature can used as diagnostic of the presence of a planet. 

2. The DyCoSS code 

In lThebaultl(l20T2h . we give a thorough description of the (then 
still unnamed) DyCoSS algorithm. Let us here only recall its 
main characteristics and the improvements that have been im- 
plemented for the planet-in-a-disc version. 

2.1. Principle 

As mentioned earlier, the basic set-up for DyCoSS is a colli- 
sionally active disc and one perturbing body. The main idea be- 
hind the code is that the grains that populate the disc at a given 
time fo, when the perturber is at a position angle <po on its orbit, 
are the superposition of grains produced at different moments in 
the past. Considering a discrete time increment Af, then grains 
present at present time to are the sum of grains that have just been 
produced (at fo); plus grains produced at to - Af, when the per- 
turber was at position angle <p-\, which have not been destroyed 
(by collisions) or ejected (by dynamical perturbations) between 
to - Af and now; plus grains produced at fo - 2Af, when the per- 
turber was at position angle <p-2, which have not been destroyed 
or ejected between fo - 2Af and now, and so on... So, if we know 
the fate of particles released at to— At, at f(>-2Af, etc., then a map 
of grain densities at fo can be reconstructed. This can be done by 
performing a series of separate runs, one with grains released 

3 an improved, more versatile version of CGA handling eccentric per- 
turbers is under development (Stark, personal communication) 
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when the perturber initally is at <po, one for 0_i , one for <f>^2, etc.. 
For each of these runs the positions of all surviving particles are 
recorded after every Af time interval, the number of remaining 
grains progressively decreasing because of dynamical ejection 
or collisional destruction. The surface density E of the disc at to 
is then obtained from 

oo 

Z(*o) = ^crfo) (f0 _ Mfl (1) 

where <x( fo)(/„-/-AO is the surface density, at to, of the run started 
at to - i ■ At. The problem is of course that this procedure would 
require an almost infinite number of separate runs, one starting at 
each initial time to-i-At. This problem can be circumvented if we 
assume that the system has reached a steady state. In this case the 
disc's profile is the same at successive passages of the perturber 
(here, the planet) at the same orbital phase. If for example, the 
planet's orbital period t or \, is equal to n or b At, then all the runs 
with the planet at initial position </>_,-, </>-,■-„„,.,,, 0-,-2«„,.,,, etc., do 
correspond to the same run, so that in practice only n or b separate 
runs are needed. 

The numerical procedure is divided into 3 steps (see TB012 
for a more complete description): 

- 1) Parent Body run. A collisionless pure N-body run is per- 
formed, where radiation pressure effects are ignored, until 
dynamical steady state is reached. The spatial distribution of 
the parent body particles is then recorded for a sample of n ol -t, 
different orbital positions of the planet on its orbit. 

- 2) Collisional runs. From each of these n or b steady-state 
parent-body discs, = 2 x 10 5 small grains are released 
following a dN oc s q ds power law. The grains' evolution is 
then followed, taking this time into account the effect of ra- 
diation pressure. Depending on its size and the local optical 
depth in the disc, each grain is assigned, at each timestep, a 
collision destruction probability that depends on grain size, 
local velocity and the local geometrical vertical optical depth 
T r fi. The relative spatial distribution of T r # is obtained from 
maps of the steady-state parent body runs and its magnitude 
is scaled by the system's average optical depth to, which is 
treated as a free input parameter. All particle positions are 
recorded at each Af = t or b/n or b interval. Runs are stopped 
once all particles have been removed by dynamical ejection 
or collision. 

- 3) Recombining. For each orbital position of the perturber, 
the dynamical + collisional steady-state of the disc is ob- 
tained by recombining the position data stored at step (2) 
following the procedure given in EquJTJ each time assum- 
ing that to is the time where the planet is at a given orbital 
position 0(Oo<i<«„„,(,- 

Results are then displayed in the form of synthetic images 
and surface brightness profiles, in scattered light, assuming grey 
scattering for all particles. They are thus valid at any wavelength 
A where the flux is dominated by sc attered light emission . Using 
the Debris Disc Simulator (DDS) of IWolf & Hillenbrand! J2005D 
for a typical debris disc located at 50 AU from its star, we find 
that the corresponding wavelength domain is A < 8/mi for a low 
mass sola r-type star (in accordan ce with the result displayed in 
Fig.7 of iKuchner & Starkl (1201 Oh ) and A < 5yum for a beta-Pic 
like A star. 

We also implicitly assume that the smallest particles in the 
system, corresponding to the radiation-pressure blowout size 
s cut , are bigger than A/In, so that all particles contribute to the 



scattered-light flux as a function of their geometrical cross sec- 
tion. Note that the absolute value for the cut-off size s cul is of no 
importance in our set-up. All simulation results can be scaled up 
or down depending on the real value of s cut , as long as the star is 
luminous enough, at least ~ 0.9M Q , so that there is such a cut-off 
size. For this smallest possible stellar mass s cut ~ ljum for com- 
pact silicates, which means that our s > A/2n condition is valid 
up to A ~ 6yum in this most constraining case, thus approximately 
corresponding to the whole domain of scattered-light dominated 
luminosity. 

2.2. Improvements 

Several upgrades have been implemented with respect to the ini- 
tial DyCoSS version. A first important update is in the way the 
collision destruction probability is estimated in the collisional 
runs. In TB012, the local geometrical optical depths r, g, which 
control the particles' collisional timescales, were derived using 
a 2-D density map that was obtained from an average of n or b 
steady-state parent body disc profiles. This map was thus in prac- 
tice azimutally averaged and only had 1-D information, an as- 
sumption that was acceptable for the external-stellar-perturber 
case where azimutal structures were limited. For an embedded 
planet, however, azimutal structures are expected to be more 
prominent and we must retain the azimutal information. As a 
consequence, we use the parent body runs to produce 100 den- 
sity maps, corresponding to 100 different orbital position^ of 
the planet, which are then used in the collisional runs (the den- 
sity map used at a given time is the one corresponding to the 
planet location closest to the present one). 

Apart from this upgrade, the collision prescription is the 
same as in TBQ12 and is gi ven by Equ. 1 of that paper. We fol- 
low [S^&^chne3 $MM) 

and assume that collisions between 
similarly-sized dust grains are fully destructive, i.e., particles are 
removed when suffering one such impact. This simplifying pre- 
scription has been taken because of numerical constraints, as a 
realistic treatment of fragmenting collisions is still impossible 
in an N-body code approach (see discussion in Section 1). It 
is of course not fully realistic, as collisions cannot fully vapor- 
ize particles on the spot and should produce clouds of smaller 
fragments. Note, however, that Equ.l of TB012 estimates the 
collision rate between similarly-sized bodies and that the typical 
impact velocities in debris discs are expected to be high enough 
for such collisions to lead to catastrophic fragmentation and thus 
produce a la rgest-remaining -fragment much smaller than the im- 
pactor (e.g., lThebaultl20 09)Flso that one can assume that a given 
impactor loses its identitys after such an impact. Furthermore, 
for the smallest grains that dominate the system's brightness in 
scattered light, of sizes b etween the radiation pressure blow-out 
size s cul and a few s cut dThebault & Augereaul 12001 . all pro- 
duced collisional fragments will be smaller than s cu , and thus 
very quickly removed from the system. As such, our prescrip- 
tion gives a satisfying estimate of a typical particle's collisional 
lifetime, which is the main parameter of interest in the present 
study. 

Another update is that Poynting-Robertson drag is now taken 
into account for the collisional runs. Finally, in order to have 
a better statistics for a wider range of grain sizes, we perform 



4 This increment is decoupled and in general much greater than n or f, 

5 This high-velocity assumption is confirmed in our runs, where the 
average collision velocities are of the order of ~ 500 m.s~' . This is much 
higher than the typical velocity, ~ 10 - 5 0m.s -1 that is r equired to catas- 
trophically destroy micron-sized grains ( Thebault 2009) 
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Table 1. Set-up for the nominal case runs. Parameters denoted 
by a * are explored as free parameters 



PLANET 

Planet/Star mass 

eccentricity 

semi-major axis 
PARENT BODY RUN 

Number of test particles 

Initial radial extent 

Initial eccentricity 
Initial surface density 
COLLISIONAL RUNS 
Average optical depth 
Number of test particles 
Size range^ 



'H = 2 x 10" 
"e p = 
a„ = 75AU 



N PB = 2x 10 5 
30 < r < 130 AU (embedded planet case 
45 < r < 75 AU (inner-disc/outer-planet 
< e < 0.01 

S oc r -1 

*T = 2 X 10 -3 
N mlm = 2 x 10 5 
i'„„ < i < 40i'„„ 



Size distribution at f = * dN(s) 



OC i' 



5 ds 



" Scut is the radiation pressure blow-out size 

* when released from the parent body population 



2 collisional runs for each given position </>, of the planet, one 
for small grains between s cut (corresponding to a ft value of 
0.5), the cut-off size imposed by radiation pressure, and 5s cut , 
and one for big grains between 5s cu , and 40s alt . These two 
runs are then merged, weighting the big-grains run by a factor 
(40s^j/ - 5s^/)/(5s^ - s 9 c ^ r ) to account for the size distribu- 
tion with slope q. 



2.3. set-up 

We consider two different set-ups. The first one is a planet em- 
bedded inside an extended debris disc, the second one is a planet 
exterior to an inner ring-like disc. We consider a fiducial "stan- 
dard" debris disc extending from 30 to 130 AU (the typical stel- 
lar distance at which most structures in resolved debris discs 
have been imaged) for the embedded planet case, and from 45 to 
75 AU for the exterior planet configuration. As for the planet, we 
take its semi-major axis to be a p = 75 AU. We consider as a ref- 
erence case a planet of mass M p /M* = jj = 2 xl0~ 3 on a circular 
orbit e p — and a system of mean optical depth to = 2 x 10~ 3 
corresponding to a dense debris disc like ft Pic. Note that the 
absolute values of the disc and planet radial positions do not 
strongly affect the results and that our results are easily scalable 
(see Sec|4]i. What matters here mostly are the disc's assumed 
geometrical optical depth tq, which controls its collisional evo- 
lution, and the planet's mass M p and eccentricity e p , which de- 
termine its efficiency in perturbing the disc and ejecting unstable 
particles. All these 3 parameters are explored as free parameters. 
All set-up parameters are summarized in TabQ] 



3. Results 

3.1. Pure N-body runs: parent body population 

In these purely dynamical runs we consider the evolution of 
collisionless large particles, not affected by radiation pressure, 
which will serve as seed particles for the collisional runs. As 
already mentioned, these collisionless simulations are only the 
first step of our complex numerical procedure and not the pur- 
pose of our study, but we will present some of them (for the 
embedded planet case only) for pedagogical purposes and also 
because they will serve as easy references to identify the effect 



of collisions and radiation pressure on the purely dynamical be- 
haviour of perturbed discs. 

For these parent body runs (hereafter PBR), we let the system 
evolve until a dynamical steady state has been reached,i.e., until 
the disc's spatial structure does not vary between two passages 
of the perturbing body at the same orbital phase (see TB012). 
For the case of a planet on a circular orbit, this happens rela- 
tively quickly, typically less than a hundred planet orbits, for all 
^ the planet masses considered (see FigQJ, b and g). For the nom- 
c |n3}l n — 2 x 10~ 3 and e p = case, several pronounced structures 
are visible. The most prominent one is the gap in the middle of 
the disc, corresponding to the chaotic unstab le zone surround- 
ing the planet's location (e.g. IWisdomlll980h . where only the 
two stable co-orbital Lagrange points L4 and L5 are populated. 
This gap corresponds to the so-called "feeding zone" surround- 
ing the planet's orbit. Clearly visible are also the 1:2 and 2:1 
mean motion resonances (hereafter MMR) at 47 and 119 AU, re- 
spectively. The apparent gaps at the resonances' location is due 
to the fact that resonant particles reach high eccentricities e and 
large radial excursions, thus spending a short fraction of their or- 
bit at the exact radial resonant location. For a lower planet mass 
fi = 2 X 10~ 4 , the gap logically becomes more narrow, and the 
whole co-orbital azimutal region is now populated, except for 
the region just surrounding the planet. 

For an eccencric planet, the system's evolution is very dif- 
ferent, and the time to reach steady state can be very long, espe- 
cially for small perturbe rs. Indeed, as has been thoroughly inves- 
tigated bv lWvattl ((2005), the dynamical response of a disc to the 
perturbing presence of a planet with e p > is the formation of 
two spiral waves that propagate outwards and inwards from the 
planet's position. As time passes, the spirals become more and 
more tightly wound so that they eventually become unnoticeable 
and appear as an asymmetric ring. The time at which this disa- 
pearance occurs depends on radial distance to the planet and is 
roughly comparable to the secular precession time (Wvatt 2005) 



tsec = 6.15 



-2.5 



1 



yb\ /2 (a p /a) j 



fjee(3:2), 



(2) 



where a p = 1 for a p < a an d a p — a p /a f or a p > a, by 2 (a p /a) is 
the Laplace coefficient (see IWvatll (120051) for more details), and 
tsecQa) is the secular timescale at the location of the 3:2 mean 
motion resonance: 



0.65 It, , 



(3) 



where t p is the orbital period of the planet. 

In agreement wifh lWvattl d2005l) 's results, we witness the de- 
velopment and progressive disappearance of these 2 spirals. The 
dynamical steady state of the system is reached when secular 
precession has rendered the spiral structures indiscernable ev- 
erywhere in the system, i.e., when t >max(f Je , c (a,„), t sec (a our )). 
For our reference case of a fi = 2 x 10~ 3 planet, and with 
a,„ = 30AUanda ouf = 130 AU this happens after ~ 2xl0 6 years. 
For a smaller /i = 2 x 10~ 4 planet, this time to steady-state in- 
creases to ~ 2 x 10 7 years (FigJTJi and f). It is important to point 
that this value exceeds the age of some bright and young debris 
discs such as HR4796 or ft Pictoris. We adopt a conservative ap- 
proach and exclude from our study all cases for which steady 
state is reached on timescales longer than this value. We show 
however, for illustrative purposes, one such case, with e p =0.1 
and /j = 10~ 5 (Fig[TJi), where the transient spirals are still clearly 
visible at ~ 2 x 10 7 years. 
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Fij>. 1. Spatial distribution at the end of the parent body dynamics-only runs, i.e., once a steady state has been reached (except for 
the yu = 10~ 5 , e\, — 0. 1 run), for different planet orbits and masses. Distances are given in AU. The planet's semi-major axis is 75 AU 
for all cases, and is aligned with the X-axis, with its perihelion on the right hand side. All snapshots are shown when the planet 
passes at aphelion (the planet's location is marked by a cross). For all runs except the ,u = 10~ 5 , e/, =0.1 one, the time given is 
the time at which steady-state has been reached. For the ,u = 1(T 5 , e h = 0.1 run, the time is the one at which simulation has been 
stopped, i.e., before steady state has been reached (see text for details). 
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Regarding the spatial distribution of the steady state disc for 
these eccentric planet cases, we see that, even for a massive ju = 
2 x 10~ 3 planet, almost no resonant structure is longer visible 
in the inner or outer discs, and the Lagrange points, although 
still present, are much less pronounced, especially in the e p = 
0.2 case. In addition, the outer disc gets strongly depleted while 
the inner disc assumes an asymmetric shape aligned with the 
planet's orbit. 

3.2. Collisional system: embedded planet 

Figf2]shows synthetic head-on images, in scattered light, at colli- 
sional and dynamical steady-state, for different planet orbits and 
collisional activities in the disc (as parameterized by its average 
optical depth to). For our nominal to = 2 X 10~ 3 case, the main 
structures observed in the pure N-body runs, such as for exam- 
ple the Lagrangian points and the 2:1 resonances in the e p = 
run, are still clearly visible, but they are much less pronounced. 
The other consequence of the coupled effect of collisions and ra- 
diation pressure is that it injects small, high-/? grains in regions 
that are devoid of parent bodies. These regions are the exter- 
nal region beyond the outer edge (at 130 AU) of the parent body 
disc, and, more importantly, the dynamically unstable region (the 
gap) close to the planet's orbit. The global effect of collisions 
and radiation pressure is thus to tend to homogenize the surface 
density profile, as is clearly illustrated for runs with higher op- 
tical depth (i.e., collisional activity), where the differences with 
the collisionless N-body simulation get much more pronounced. 
For the e p = orbit, the accumulation around the Lagrangian 
points is for example replaced by a horseshoe co-orbital ring for 
higher value of tq (Fig[2}i). Conversely, for lower values of to, 
the spatial profile becomes much closer to that of the parent bod- 
ies (Figf2g). 

Despite this blurring effect of collisions, for all the head-on 
images displayed in Figf2] the planet's radial location remains 
relatively easy to determine from the position of the radial gap. 
For high-mass low e p cases, even the planet's azimutal position 
remains easy to spot in between the well defined Lagrangian 
points, with the noticeable exception of the high to case (Figf2ji) 
where the steady production of small collision fragments erases 
almost all azimutal structures in the co-orbital region. 

The situation becomes very different if the system is seen 
edge-on. Figj3] shows the midplane luminosity profile along a 
radial cut passing by the planet's current position. While the dis- 
tribution of large parent bodies shows a pronounced gap around 
the planet's radial position (FigO, with the flux dropping by a 
factor ~ 5 in the 50-90AU region, the coupled effect of colli- 
sions and radiation pressure on small grains greatly alters this 
picture. Even for a relatively low to = 4 x 10~ 4 value, the lu- 
minosity drop around the planet falls to less than a factor 2. It 
is of less than 30% in the nominal To = 2 x 10~ 3 case, and falls 
down to a few percent for even more collisionally active discs 
(to = 0.01). Furthermore, this small flux variation is not peaked 
at the planet's location but diluted over a rather extended region, 
so that the planet's effect on the profile is a wide and poorly de- 
fined knee. 

The effect of the planet on the luminosity profile is even 
weaker when considering masses smaller than the fi = 2 x 10~ 3 
super- Jupiter of our reference case. As can be clearly seen in 
FigS for a dense debris disc (to = 2 x 10~ 3 ), a Saturn-like planet 
(ji — 2x 10~ 4 ) leaves a barely perceptible signature, and a super- 
Earth (ji = 10~ 5 ) results in a radial profile almost identical to 
the no-planet case. This is in stark contrast to the profile that 
would have been obtained had collisions and radiation pressure 
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Fig. 3. Synthetic midplane luminosity profiles for a disc seen 
edge-on, in scattered light, for the standard pi — 2 x 10~ 3 and 
e p — case and for different values of the disc's average optical 
depth (i.e., collisional activity). The shown profiles correspond 
to a cut along the X axis when the planet passes at its maximum 
elongation along this axis (the planet's location is marked by a 
black circle). Profiles are computed assuming a grey scattering 
function. The profile for the large parent bodies is directly de- 
rived from the PB runs presented in Section [3~Tl 
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Fig. 4. Same as Fig|3] but with e p = and t = 2 x 10~ 3 and 
when varying the mass of the perturbing planet. 



been neglected. Fig|5]shows such profiles, derived from the par- 
ent body runs. As can be seen, in these fiducial cases, there is 
always a well defined and narrow gap close to the planet's radial 
location, even for the smallest case of a fj. = 10~ 5 super-Earth. 
This clearly illustrates how the coupling of collisional activity 
and radiation pressure efficiently erases the signature of a planet 
(see Sec|4]for further discussion). 
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urations (see text for details). 
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Fig. 6. Same as Figj3] but for different eccentricities of the 
planet's orbit (its mass being fi = 2 X 10" 3 ). 



Fig. 7. Inner-ring/outer-planet configuration. Steady-state syn- 
thetic surface brightness, when the system is seen head on, along 
a radial cut passing by the planet's current location. The grey 
area marks the initial radial extent of the parent body ring. 
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The blurring effect of collisions is even more apparent for 
an eccentric embedded planet. For the e p =0.1 case, collisional 
production of small high-/? grains homogenizes the region co- 
orbital with the planet (Figf2^), in sharp contrast to the pure 
N-body runs where the co-orbital region is almost empty, safe 
for the 2 Lagrangian points. The same effect is observed for the 
e p = 0.2 case, where no azimutal inhomogeneity is visible in 
the co-orbital region on the 2-D map (Figf2}l). For this specific 
case, however, the collisional runs confirm the main result of the 
parent-body runs, i.e., a strong depletion in the region beyond 
the planet's orbit. This depletion is so strong that it almost com- 
pletely erases the frontier between the planet's co-orbital region 
and the outer disc. If the system is seen edge-on, then the sig- 
nature of an eccentric planet becomes even weaker (Fig|6|l. It is 
even less pronounced than in the e p = case. For the e p = 0.2 



Fig. 8. Inner-ring/outer-planet configuration. Relative contribu- 
tion of the "big" grains (J3 < 0.12) to the total luminosity for 3 
different cases. 



run, there is almost no variation in the radial luminosity profile at 
the planet's distance. Here again, this result is in sharp contrast 
to the distribution of large parent bodies, for which, regardless 
of the value of e p , a clear dip is always visible in the edge-on 
luminosity profile. 

3.3. Collisional system: Inner disc/external planet 

We now consider a different configuration, where the planet is 
located exterior to an inner debris ring. The main objective is 
here to see to what extent the planet can shape the outer edge of 
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Fig. 9. Same as Fig[7J but for different eccentricities of the exte- 
rior planet (and for fi = 0.002 and to = 0.002). For the e p > 
cases, plots are presented for a passage of the planet at periastron 
and apoastron. Note that, since for these e p > cases the unsta- 
ble chaotic region inside the planet's orbit is much wider than 
when e p = 0, the initial disc of parent bodies extended down to 
25 AU for the e p = 0.1 and e p = 0.2 runs. 

the disc, and in particular if it can affect the luminosity profile in 
oc r~ that should "naturall y" form beyond the radial location 
of a collisional debris ring (Theba ult & Wull2008h . We consider 
a standard set-up where the inner ring of parent bodies initially 
extends from 45 to 75 AU and the planet is positioned at the outer 
edge of the ring, i.e., 75 AU. We restrict ourselves to the case of 
a planet on a circular orbit, and explore its mass and the disc's 
collisional activity as free parameters. 

Fig|7] shows a radial cut of the surface brightness profiles 
when the system is seen head-on. The first effect of the planet 
is to truncate the disc of parent bodies, at a distance corre- 
sponding to its chaotic zone. For each different curve, the lo- 
cation of the outer edge of the stable parent body ring corre- 
sponds approximately to the transition point where the bright- 
ness peaks before rapidly decreasing with radial distance, i.e., 
at r ~ 61 AU for a fi = 2 x 10~ 4 planet, ~ 52 AU for a 
fi = 2 x 10~ 3 planet and ~ 50 AU for a fi = 4 x 10~ 3 one. 
Interestingly, beyond this outer edge rpg of the PB disc, the 
brightness profile has a slope that only weakly depends on the 
planet's mass. For the no-planet run, beyond a transition re- 
gion extending from 75 to ~ 95 AU, it assumes the standard - 
3.5 value due to the steady collisio nal production of radia tion- 
pressure sensitive small grains (see iThebault & Wull2008l) . For 
a low-mass planet of fi — 2 x 10~ 4 this profile is basically un- 
changed, and it only gets slighly steeper, ~ -4 to -4.2 for more 
massive planets (pi = 2 x 10~ 3 and fi — 4 x 10~ 3 ). 

These results might seem counter-intuitive, because the r > 
rpB domain is by definition dynamically unstable and we would 
expect the planet to eject a large fraction of the small grains, 
produced in the main ring, that populated this outer region. 
However, this ejection process takes time, basically needing a 
close-encounter with the planet, so that small grains can have a 
long time of residency in this dynamically "forbidden" region. 
True, this survival time is shorter than for the no-planet case and 



we would expect to see less grains beyond rpg than in this un- 
perturbed case. But this depletion of very small fragments in the 
outer regions is partly compensated by the injection of slighly 
bigger collisional fragments due to planetary perturbations. This 
is clearly illustrated in Fig[8] showing that /3 < 0. 12 grains (with 
size s > As cut ) produced in the PB ring, which have an apoas- 
tron < rpslil - 2/3), and should in principle not be present in 
the r > rpslil - 2/3) region, contribute to up to ~ 20% of the 
brightness in these regions in the presence of a Jovian perturb- 
ing planet. 

Profiles close to the reference -3.5 value are also obtained 
when considering planets on eccentric orbits (FigO. In these 
cases, the stable ring of parent bodies is truncated much further 
inside the planet's orbit (for the e p = 0.2 case, it only extends 
to ~ 3 1 AU and ~ 37 AU in the periastron and apoastron di- 
rections, respectively) but the profiles beyond this main ring are 
comprised between 3.4 and 3.9. 

The only way to obtain a significantly steeper brightness pro- 
file is to assume a much lower collisional activity within the disc 
(to = 2 x 10~ 4 ). In this case, the balance between collisional 
injection of small grains and dynamical ejection by the planet 
is shifted in favour of the latter, and the slope of the luminosity 
profile is 5. 

When considering head-on images of the disc (FigfTUt. it can 
be seen that, despite of the blurring effect of collisions coupled 
to radiation pressure, the planet's perturbations can induce sig- 
nificant spatial structures, mostly in the region of the parent body 
ring. For a fi = 0.002 perturber, the width of this main ring dis- 
plays azimutal inhomogeneities, with three wider "blobs", lo- 
cated at ~ +30 and 180 degres of the planet's location, that pre- 
cess with the planet. Not surprisingly, these strucutures are more 
prominent for the low collisional activity case (to = 0.0002), 
but are still visible even in the nominal to = 0.002 case. For an 
eccentric planet, these azimutal structures disapear, but the main 
ring takes on an eccentric shape aligned with the planet orbit's 
main axis. From the location of the peaks on the apoastron and 
periastron sides of the ring's radial profile (Fig|9j, we estimate 
the ring's eccentricity to be e r i ng ~ 0.08 for e p = 0.2. This value 
is consistent with the forced eccentricity that can be derived from 
the secular theory of Laplace-Lagrang e for this u value and this 
disc-to-planet relative distance (e.g. iMurrav & Dermott| [l999: 
iMustill & Wvattl2009l) . Another clearly visible feature is the fact 
that the periastron side of the ring is brighter than the apoastron 
one. T his is the pericentre glow effect identified by I Wvatt et ail 
( 1999), due to the fact that the periastron side of the ring, being 
closer to the star, receives and scatters more light. From Fig|9] 
we see that the pericentre side is ~ 1.35 times brighter than the 
opposite one for the e p = 0.2 case, which compares well with the 
predicted difference in received light, i.e., ((1 +e n „g)/(l -e r i„ g )) 2 
for the measured e rmg = 0.08 value. 

For a lower mass planet (fj. = 0.0002), no significant spatial 
planet-induced structures are visible on the head-on images. 

4. Discussion 

The main question the present study aims to address is if planets 
can significantly affect debris disc structures and more specifi- 
cally if the presence of planets can be indirectly inferred from the 
spatial structure of an observed disc. As noted in Section 2, we 
here implicitly consider observations in scattered light where the 
luminosity is dominated by the smallest grains in the disc, i.e., 
for wavelengths shortward of ~ 5 - 8 //m. This is, however, not a 
too stringent limitation, since a large fraction of resolved images 
of debris discs, and in particular those with the best resolution, 
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Fig. 10. Inner ring/outer planet: head-on synthetic images in scattered light, at steady state, for four different set ups. The graphs 
on the left hand side show the disc profile when the planet passes at apoastron, and the graphs on the right hand side corre- 
spond to periastron passages. The dotted curve shows the planetary orbit. An animated version of these graphs can be found at 
http://lesia.obspm.fr/perso/philippe-thebault/planpert.html 
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have been obtained in scattered light. The situation might, how- 
ever, change wit h the upcoming h igh angular images expected 
from ALMA (see lErtel etai]|2012l) . 

Another implicit assumption is that there is a minimum size 
s cut below which small grains are blown out by radiation pres- 
sure. This implies that the central star has a mass of at least 
~ O.9M . Note, however, that for lower mass stars of type M, 
stellar wind could play a role similar to that of radiation pres- 
sure around massive stars, blowing away grains below a given 
threshold size. This is for instance what is expec ted to take place 
in the AU Mic system (lAueereau & BeusfeOOrJ) . 

Considering the already huge parameter space for the prob- 
lem, all simulations were carried out for one given semi-major 
axis of the planet a p = 75 AU. However, our results can be 
easily scaled to other planet separations with the same relative 
disc/planet configuration, provided that the system's optical 
depth to is the same. In this case, the collisional timescales fol- 
low the same scaling laws as the dynamical ones (those linked to 
the planet), i.e., r ■ . The only difference between our reference 
case, with a p = 75 AU and a 30 < r < 130 AU disc, and, say, 
a system with a a p — 25 AU planet and a 10 < r < 43 AU disc 
is that all timescales should be shortened by a factor ~ 5. Note, 
however, that there is a limit to arbitrarily scaling up our results 
to high values of a p , which is that impact velocities in the disc 
should remain high enough to lead to erosive collisions. Apart 
from this global timescale effect, the three main parameters that 
control the system's evolution are the planet-to-star mass ratio fi, 
its eccentricity e p and the disc's collisional activity, as controled 
by its average vertical optical depth to. All these parameters have 
been explored in our numerical investigation. 

4. 1 . Dectability of an embedded planet in scattered light 
images 

4.1 .1 . Head-on viewed disc 

Our simulations have shown that the coupling of collisional ac- 
tivity and radiation pressure can efficiently attenuate or even 
erase some spatial structures that an embedded planet might cre- 
ate within a debris disc. This is particularly true in the feed- 
ing zone surrounding the planet's orbit. This region appears as 
a deep and well defined gap in the distribution of the parent 
body population, with the exception of the two Trojan popu- 
lations at the Lagrangian points, but takes on a more homo- 
geneous aspect once collisional production of small radiation- 
pressure-affected grains is turned on. This global tendency of 
co llisions to attenuate sp atial structures has also been witnessed 
by [Stark & Kuchnerl d2009l) . but for a different set-up, where the 
planet lies interior to a ring that migrates inwards because of 
Poynting-Robertson drag, and for smaller planet masses. For the 
present case of an embedded planet, and even with massive, 
super-Jovian planets, azimutal structures like the Lagrangian 
points are no longer visible and the feeding zone appears as an 
almost continuous ring. However, the radial structure of the disc 
still bears the scars due to the presence of the planet: the feed- 
ing zone surrounding the planet's orbits is, for most cases, still 
less densely populated than the rest of the disc and appears as 
a gap-like structure. As a consequence, on 2-D face-on images, 
the presence of the planet, and even its radial location, can eas- 
ily be inferred from this ring-gap-ring structure, and this even for 
Saturn-like or super-Earth planets. Although, in the latter case, 
the width of the gap, i.e., ~ O.Q5a p (Fig|2};), might be too small 
to be observationally detectable. In addition, for the specific case 
of a Jovian planet on a circular orbit (Fig|2^), the structures due 



to the 2:1 resonance are still visible enough in order to also in- 
fer the azimutal position of the planet. The only cases for which 
the planet does not leave a clearly identifiable gap in a head-on 
viewed collisional disc is if the planet has a high mass and is on 
an eccentric orbit. In this case, the region exterior to the planet's 
location is strongly depleted and its contrast relative to the feed- 
ing zone is very low. However, in this case the planet's influence 
might be inferred from another feature, i.e., the eccentric shape 
of the inner disc (Fig|2}i). 

4.1.2. Edge-on viewed disc 

For a system seen edge-on, the situation is much less favourable. 
Figsj3]|4]and|6]clearly show that the signature left by the planet 
in the disc's radial luminosity profile is in most cases very weak. 

We firstly note that there is almost no asymmetry between 
the two opposite sides of the disc, i.e., the luminosity profile of 
the side where the planet is located is almost indistinguishable 
from the profile of the opposite side. Even for a perturber on a 
moderately eccentric orbit (e p < 0.2), the differences between 
the two ansae remain relatively limited, in any case too limited 
to be used as a smoking gun proof of the planet's presence. 

Secondly, contrary to the parent body distribution, where 
there is a narrow and pronounced luminosity drop at the planet's 
location even for low mass planets (Fig|5]l, in a collisional disc 
with small grains this sharp drop is replaced by a much wider, 
and often poorly defined and shallow hollow. For most explored 
parameters, the main effect of the planet is not this shallow hol- 
low but a steepening of the luminosity profile in the radial re- 
gion interior to the planet's position, as is best illustrated in 
the case of a /i = 2 x 10~ 3 planet on a e p = 0.2 orbit (Fig|6]l. 
Given the spatial resolution obtaine d for the best resolved ed ge- 
on debris disc s, such as ffPic (e.g ., iGolimovski et al.ll2006l) or 
AUMic (e.g., Met chev et al.ll2.Q05h . this steepening should po- 
tentially be observable, at least for bright and nearby systems. 
However, even if it is observed, it is doubtful that such a steeper 
profile could be unambiguously associated to the presence of a 
planet. Indeed, if on Fig|4]the difference with the planet-less case 
is clearly visible, for real systems there would of course be no 
such "reference" case showing what the system should look like 
without a planet. Without this reference case, a steep slope could 
as well be associated to an intrinsic steeper radial distribution of 
the parent bodies (planetesimals) that could mimic the potential 
effect of a planet. 

In this respect, the small hollow that the planet induces in 
the profile would be a less ambiguous observable signature, pro- 
vided it is pronounced enough to be observed, because it would 
be harder to mimic in a planetless disc where collisional pro- 
duction of small grains should naturally even out any short-scale 
variations in the density (i.e., luminosity) profile. If we consider 
that a hollow is "pronounced" when there is an inflection point in 
the luminosity profile, then we see that this criteria is just about 
met for our nominal /j — 2 x 10~ 3 , e p = and to = 2 x 10~ 3 
case (Fig f3j. Smaller planets (Saturn-mass and below) would 
induce much shallower hollows, as would also planets (even 
massive ones) on e p > 0.1 orbits. The same is true for denser 
discs (to = 0.01) where the intense collisional activity elimi- 
nates the inflection in the radial luminosity profile. Conversely, 
the hollow is logically more pronounced for more tenuous discs 
(to = 4 x 10~ 4 ). However, it remains to see if, even for the lim- 
ited parameter space where the hollow is "pronounced", it would 
be concretely observable in real discs. For one of the best re- 
solved edge-on discs, i.e., the one around AU Mic, such narrow 
dips seems indeed to be observed twice in the NW arm, once 
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aroun d 25 AU and once at ~ 50 AU (see Fig.4 of iMetchev et al.1 
I2005h . with both times a depth that exceeds the plotted error bars. 
However, the presented error bars correspond to lcr only. For 
more reliable 3cr error bars, none of the dips would remain vis- 
ible. Furthermore, in both cases, these dips have no counterpart 
in the opposite SE arm, whereas our simulations predict that a 
planet-induced dip should be observable on both sides. Thus, we 
rate these structures as noise and conclude that the predicted in- 
flection point is probably not observable directly with present 
instruments. 

On the other hand, our simula tions show that a cle ar jump 
in the profile (not just a knee as in lMetchev et al.ll2005l) . with a 
well defined, similar slope inside and outside the inflection point 
could reveal the presence of the inflection point without the need 
of detecting it directly (a local flattening of the profile as found 
for our to = 0.01 case in Figj3]is already sufficient to produce the 
effect). Such a feature might be observable, because the slope of 
the surface brightness over a broad range of radii is usually well 
constrained as it is sampled by a large number of data points. 
Nevertheless, it is important to note that this is only a qualitative 
statement on the example of one observation. The exact sensitiv- 
ity depends on a large number of parameters of the observation 
(e.g., integration time, PSF subtraction technique) and the disc 
(e.g., surface brightness, extent) and has to be evaluated from 
case to case. A full, quantitative discussion is beyond the scope 
of this work. W e note, however, tha t near-future instruments like 
VZjy SPHERE dDohlen et al.ll2006l) are expected to significantly 
exceed present instruments in terms of sensitivity and contrast 
(i.e., the ability to observe fainter discs and/or at lower angular 
separation from the star). This will allow to obtain data of similar 
or increased quality of a larger sample of debris discs, improving 
the chances to detect the predicted structures. 

Regardless of these considerations, it must be pointed out 
that, even if an inflection point is observable in the luminosity 
profile and if it can be linked to the presence of a planet, it does 
not directly give the radial position of the planet, since it is al- 
ways located further out from the star. The offset between the 
planet's radial location and the position of the inflection point 
depends on the planet's mass and on the collisional activity in the 
disc. It varies between ~ 0.05a p and ~ 0.2a ; , for the parameters 
explored in our simulations (see Figs|3] and Unfortunately, 
the value of the offset cannot be used to infer the planet's char- 
acteristics because the problem is degenerated: there is a trade- 
off between the planet's mass and the disc's collisional activity, 
which can both affect the amplitude of the offset. 

4.2. Ring sculpting by an outer planet 

4.2.1 . Can a planet truncate an inner ring? 

The dynamical sculpting of an exterior perturber has been of- 
ten invoked as being a possible explanation for debris discs, 
or debris rings, with sharp outer edges, "sharp" meaning here 
sharper than the "natural" luminosity slope in r~ 35 beyond the 
main ring's outer li mit r o ut (The bault &Wull2008h . 

iThebault etaD (l2010h and TB012 have shown that a stellar 
companion can to some extent significantly deplete the region 
beyond r oul , even if it is never able to fully truncate a collision- 
ally active disc. The most effective depletion is obtained for a 
binary on a circular orbit lying as close as possible to the ring, 
i.e., so that the outer edge of the parent body ring corresponds to 
the limit beyond which orbits become unstable (in other words, 
the binary truncates the parent body ring at r = r out ). In this 
case, the flux at a distance 2r our from the primary is decreased 



by a factor 10 (with respect to the flux in the no-planet case) in 
the presence of a fx = 0.5 companion star on a circular orbit, 
and the slope can be as steep as ~ r~ 8 in the [r out , 1.5r out ] re- 
gion, before slowly decreasing again towards 5 at larger dis- 
tances. For more eccentric binaries, the depletion effect is more 
limited, mainly because the companion spends much time at its 
apoastron far from the ring, thus leaving time for collisions in 
th e ring to populate the "forbidden" r > r out domain (see Fig. 2 
of IThebault etaDIIoTob . 

The present simulations show that planets are less effective 
than companion stars for truncating collisionally active debris 
rings. For a dense to = 2 x 10 -3 ring, the slope of the luminosity 
profile beyond r out is relatively close to the -3.5 standard value in 
the planet-less case, and this regardless of the perturbing planet's 
mass. A maximum slope of only ~ -4.2 is reached for a. p - 
4 x 10~ 3 perturber (Fig|3. 

Note, however, that these slopes are not reached immedi- 
ately outside the parent body ring, but only after a small tran- 
sition region of relative width Ar/r ou , ~ 0.2, where the profile is 
steeper. This is in particular true in the planet-less case, where 
the slope of the luminosity profile is ~ -7 in the [r out , 1.2r OM ] 
region. As a matter of fact, for small planets (p < 2 x 10~ 4 ), the 
profile in this narrow region is less steep than in the no-planet 
case, mainly because planetary perturbations inject medium- 
sized grains (J3 ~ 0.1 - 0.2) in these regions just beyond r ou , 
(Fig|8]l. For more massive planets, the profile in this transition 
region gets steeper, but not much steeper than in the planet-less 
case. This point is of importance for the specific case of the 
HR4796 ring, for which preliminary DyCoSS r esults were pub- 
lished in a recent paper dLagrange et al.ll20 1 2h . The main con- 
clusion of that paper, i.e., that the steep luminosity profile in 
~ r~ 9 outside the main ring could only be fitted with a massive 
outer planet has to be reinvestigated. Indeed, the region where 
the S/N of the observed profile is high enough does not extend 
very far outside the main ring (see Fig. 16 of that paper), i.e., 
not far outside the narrow transition region where slopes could 
be high even for a planet-less case, a case that was not consid- 
ered in Lagrange e t al.l d2012l) . We present in Fig[TT] a revised 
version of that paper's Fig. 16, taking this time into account the 
no-planet case. As can be seen, in the 75-90 AU region where 
the observed profile has a reliable S/N ratio, the no-planet case 
does in fact give a better fit than the "small" planet case (3Mj up ). 
However, the massive planet (8Mj up ) cass is still the one giving 
the best fit, so that the main conclusion of Lagrange et al. (2012) 
still holds. But to unambiguously discriminate between the dif- 
ferent scenarios, reliable observations further outside the main 
ring are clearly necessary. 

As a matter of fact, there is to our knowledge no disc for 
which a steep (< r~ 5 ) luminosity profile has been observed far 
outside the expected location of the main parent body ringQ. For 
the two archetypal examples of rings with steep luminosity de- 
creases, Fomalhaut and HR4796, these profiles have not been 



6 Note that, since HR4796A is a massive A0V star of mass 2.2 M Q , 
the 8 and 3Mj„ p cases correspond to p = 0.0035 and p = 0.0013, re- 
spectively 

7 We here define the outer limit of the main ring as the outer limit 
of the large parent bodies ring, which roughly corresponds to the point 
where the scattered light luminosity peaks. This outer limit is thus dif- 
ferent from the one usually given in observational studies, i.e., the radial 
distance where the flux falls to half that of the peak. From our study we 
find that this observational outer limit is actually located beyond the 
outer edge of the main collisionally active ring. For most cases, this 
limit is close to the outer edge of the transition region, with steeper 
slopes, just outside the main PB ring 
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Fig. 11. Head-on viewed surface brightness profiles for the spe- 
cific case of the HR4796A disc. The parent body ring is located 
between 67 and 75 AU. For each case, the planet is placed as 
close as possible to the parent body ring (i.e., so that the outer 
limit of the stability region coincides with the ring's outer edge). 
The s olid line indicates the d eprojected observed profile, derived 
from ISchneider et all J2009). The horizontal line delineates ap- 
proximately the part (abo ve this line) of the pro file with a dy- 
namical range of 10. As in Lagrange et al. (2012), we show, for 
each case, the radial cut at the "best" position angle, i.e., the one 
providing the best match to the observed profile. Note that the 
goal of these runs is to fit the outer profile, not the inner one, 
which sho uld probably be shape d by other mechanisms (see dis- 
cussion in Lagrang e et al.ll2012l) . so that not too much physical 
significance should be given to the profile interior to the ring's 
inner edge. 



reliably obtained beyond the expected transition region just out- 
side the main ring. 

From our simulations, the only way to get profiles that are 
significantly steeper than in the planet-less case, both in the nar- 
row transition region outside the PB ring and further out, is when 
assuming a much lower collisional activity in the disc. This is 
clearly illustrated by the To = 2 X 10~ 4 run in Figj7] for which 
the slope in the transition region is ~ -10 and tends towards 
~ -5 further out. Such a low ~ 2 x 10~ 4 optical depth is, how- 
ever, much lower than the one derived for HR4796, ~ 5 x 10~ 3 
dAugereau et all 119991: IWvatt et all 1 19991) . Nevertheless, t ~ 
2 x 10~ 4 values are comparable to those de rived for other re- 
solved debris ri ngs, in particular Fomalhaut dBolev et al.ll2012l) 
and HD202628 dKrist et alj2012l) . The ~ -5 value we find could 
be, given the uncertainties on the data, compatible with the mea- 
sured slopes in the outer regi ons of these 2 sys tems, i.e., ~ -4.6 

to 6.5 for FomalhauFT dKalas et al.ll2005l) and 4.7 for 

HD202628. However, the radial extent of the reliable data for 
Fomalhaut does not exceed ~ 1.15r„„, and is thus confined to 
the "transition" region mentioned before, so that no conclusions 
can be drawn for this system. For HD202628, on the contrary, it 
reaches far outside this transition region, i.e., ~ l.5r out , so that 

8 The value given in Kalas et al. (2003) is -4.6, however, from the 
data given in Fig. 3 of that paper we rind a best fit value of 6.5 



this system could in principle be explained by the presence of 
an outer jovian planet. However, given that the radial luminosity 
of this systems peaks at ~ 190 AU, the planet would have to be 
located at ~ 280 AU (as implied by rescaling Fig|7]), which is 
much further out than any exoplanet detected to date. 

4.2.2. Ring morphology 

Another important issue is the spatial structure of the main par- 
ent body ring itself. Our simulations show that massive planets 
on circular orbits can induce azimutal inhomogeneities that pre- 
cess with the planet. These inhomogeneities remain visible in 
scattered light even for a high collisional activity within the ring, 
albeit with lower contrast (FigfTUt. 

As of today, no such intrinsic brightness inhomogeneities 
have been unambiguously observed in the 7 ring-like discs 
that have been imaged in scattered light (see Section 4.6 of 
Krist et al.ll2012l) . apart from those that might depend on scat- 
tering phase angle. Some bright spots have been observed in the 
HR4796A ring, but due to the limited sig nal to noise of the data , 
it is possible that it could be an artifact (Thalm ann et al.ll201 lb 
Lagr ange et al.ll20l2l) . Note also that the detection of such pre- 
cessing inhomogeneities would not be a smoking-gun proof of 
the presence of an outer planet, as an inner planet would also 
produce such precessing structures. There is, however, a way 
to discriminate between these two hypothesis, by measuring the 
precession rate of the structures: if it is faster than the orbital pe- 
riod at the ring's distance, then it is due to an inner planet, if it 
is slower then it is due to an outer perturber. Depending on the 
quality of the precession measure, the orbital period and thus the 
semi-major axis of the planet can in principle be estimated. 

For planets on eccentric orbits, no precessing azimutal struc- 
tures are visible in our results, but the main ring becomes itself 
eccentric, with e rmg close to the secular forced eccentricity. We 
observe the well known pericentre glow effect, i.e., a brightness 
asymmetry between the periastron and apoastron sides of the 
ring. Our simulations show that this effect remains observable, 
in scattered light, even in collisionally active systems. The for- 
mation of an eccentric ring in response to an outer planet's per- 
turbations is an interesting outcome, as several such eccentric 
rings have been observed, in particular around Fomalhaut and 
HD202628 dKrist et alj|2012l) . We stress, however, that the pres- 
ence of an outer planet alone cannot explain all the characteris- 
tics of these rings. This is particularly true for the inner edge of 
these rings, which is always very sharp and has to b e sculpted by 
"som ething", most probably an inner planet (e.g. IChiang et al.l 
2009). If this inner planet is on an eccentric orbit, then the ring 
naturally assumes an eccentric shape without the need of an ad- 
ditional external planet. The question is then if there are other 
ring characteristics that cannot be explained without invoking 
such a planet. Of course, an outer perturber could easily explain 
the sharp outer edge of the parent body ring. Such a sharp edge is 
indeed a prerequisite to obtain a scattered light brightness profile 
decreasing as r~ 3 5 , or even st eeper, beyond it. F or Fomalhaut, re- 
cent 850 fim ALMA images ( Bole v et al.l2012l) provide unprece- 
dented information about the spatial distribution of the large par- 
ent bodies and sho w that they indeed form a very confined ring. 
IBolev et aD, d2012l) argue that such a confinement can be pro- 
duced by two Earth-to-super-Earth shepperding planets. There is 
however a potential problem with this hypothesis when consid- 
ering the brightness profile in scattered light. Indeed, our sim- 
ulations show that, even for the low to ~ 10~ 4 optical depth 
estimated for Fomalhaut, an Earth-mass planet cannot prevent 
the brightness profile from being ~ -3.5 beyond the main ring. 
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This seems to be in contradiction with the steep slopes, between 
-4.6 and -6.5, that have been observed in scattered light images. 
However, let us again stress that these profiles have not been 
reliably obtained beyond ~ 160 AU, less than 20 AU beyond the 
main ring's outer limit at ~ 140 AU, which is still inside the tran- 
sition ~ [r out , l.2r OM ] region where steep slopes can be obtained 
regarless of the presence of a planet (Fig|7]l. Additional scattered 
light observations, probing the Fomalhaut disc further out from 
the main ring, are needed in order to confirm or invalidate the 
outer-planet hypothesis. 

An alternative explanation for a sharp outer edge of the par- 
ent body ring could be that this is the natural outer limit beyond 
which planetesimals could not form during the proto-planetary 
disc phase. Our simulations show that it would be difficult to dis- 
crimate between such a system, with an inner eccentric planet 
close to the "natural" outer limit of a debris disc, and a system 
with the same inner planet and an outer planet on an eccentric 
orbit truncating the disc. We note that an even more radical alter- 
native scenario, dispensing even with the inner planet, has been 
recently proposed to explain the formati on of sharp and eccentri c 
debris rings in a planetless environment dLvra & Kuch ner 2012). 
This issue is thus still an open one. 

4.3. Comparison with previous works and perspectives 

As mentioned in Sec.l, the present study is by no means the 
first numerical investigation of the signature left by a planet 
in a debris disc. It is, however, one of the first to take into 
account the crucial role played by collisions within the disc, 
whereas most past studies were based on collisionless dynam- 
ical s i mulations of size- less particles (e.g., iKuchner & Holmanl 
12003b iReche etal.l l2008). or, in more advanced versions, a com- 
bination of such collisionless simulations for single-sized par- 
ticles that were later recombined assuming a fiducial "col- 
lisional equilibrium" size distribution for all disc par ticles 
(lOuillen & Thomdikel2002t IWvattll2006b lErtel et aUl2012l) . 

As such, it is not surprising that our results sometimes sig- 
nificantly depart from those obtained in pure N-body studies. 
As an example, we compare our res ults to thos e obtained by 
the most recent of such studies, that of Ertel et al .1 (120121) . which 
produced synthetic images of planet-bearing discs from visible 
to far-IR and millimetre wavelengths. As mentioned earlier, at 
such long wavelengths the flux is no longer dominated by the 
smallest grains and our modelling is not valid. Thus, we only 
compare our results with their ones up to A — 10/im. One im- 
portant conclusion of lErtel etal.1 (120121) was that Jovian planets 
always induced prominent structures in discs in scattered light, 
and that such structures might in most cases be observable with 
near-future facilities, at least in face-on oriented discs. When 
w e compare our di sc structures (Figj2j> with the ones shown 
in lErtel et aTl d2012l) for the face-on case, we find that the res- 
onant structures found there to appear due to transport mech- 
anisms (Poynting-Robertson drag and stellar wind drag) dis- 
appear when considering collisions. In contrast, we observe a 
horseshoe structure and a c lear gap in the disc, which were ob- 
served by Erte l et all (1201 2l) . but only at longer wavelengths and 
not in scatt ered light. This is because the absence of collisional 
lifetimes in lErtel et al.l d2012l) allowed Poynting-Robertson drift 
to populate the gap around the planet's location even with large 
grains. We do not observe this inward drift in the present simula- 
tions, even for the lowest collisional activity considered, because 
the collision destruction rate of large grains is always higher than 
typical timescale for PR drag. As a matter of fact, the distribution 
of the biggest grains considered in our runs, of size s = 40s cut , 



is undistinguishable from that of the parent body simulations. In 
the present runs, the grains that are present in the resonant and 
horse-shoe structures were not captured when drifting inward 
from an outer region, but were present in these regions from the 
beginning. 

lErtel et ail (1201 2|) also obtained spatial structures for edge- 
on seen discs, albeit less pronounced than for head-on systems. 
Our simulations have shown that such structures can probably 
not remain in a collisionally active system. 

The only other debris disc model that incorporates a degree 
of collisional processes into an N-body scheme is the CGA of 
IStark & Kuchn er (2009). It is, however, difficult to compare re- 
sults obtained with both codes since they have not been designed 
to address the same perturber/disc configurations and have so far 
been used for different set-ups. CGA has been used to investigate 
the case of a perturbing planet interior to an outer ring, where 
an important effect is the slow drift due to PR-drag from ring 
particles into the planetary region. This drift results in efficient 
trapping into outer mean motion resonances and thus azimutal 
structures much more pronounced than for non-drifting parti- 
clefl On the contrary, PR drag has a negligible effect for the 
two set-ups that have been considered in the present study, i.e., 
embedded planet and inner-ring/outer-planet. In the former, col- 
lisional timescales within the disc are always much shorter than 
PR drag ones. In the latter, PR drag is observable, but its main ef- 
fect is to progressively move small grains into the region interior 
to the ring's inner edge (see for instance FigfTTI). i.e., far from the 
region of interest between the planet and the ring's outer edge. 
Note, however, that our results for the embedded-planet case re- 
garding the insignificance of inward d rift of big grains t oward s 
the planet is consistent with that of Kuc hner & Starkl (|2010), 
who did not observe any drift of large particles from an outer 
belt towards their inner planet. In both cases, finite collisional 
lifetimes do pr event this drift from happ ening. The apparent dis- 
crepancy with IKuchner & Starkl d2010i) 's synthetic images (see 
their Fig. 9), i.e., the presence of big grains trapped into horse- 
shoe structures and MMRs, simply follows from different initial 
set-ups: a disc-embedded planet in the present runs, as opposed 
to a planet lying in an empty region well interior to an exterior 
belt. In addition to this difference in the planet-to-disc configu- 
rations, DyCoSS and CGA have also considered different plane t 
masses. IStark & Kuchnerl (120091) and IKuchner & Starkl OOlOl) 
have explored perturber masses in the Earth-to-Neptune mass 
range, while TB012 and the present study mostly focused on 
stellar perturbers and massive, Saturn to super- Jovian planets. 

Finally, we want to underline that DyCoSS (as CGA) is not 
a fully integrated dynamics+collisions model. Although some 
important issues can be investigated with it, in particular those 
related to timescales and production rates, the treatment of col- 
lisions relies on simplifying assumptions and is not fully incor- 
porated into the scheme of the code. Furthermore, only systems 
that have reached dynamical and collisional steady state can be 
studied. It is not suited to investigate transient or violent events 
that might take place in many debris discs. As emphasized in the 
Introduction, fully-integrated codes are still out of reach today 
but improvement in computational capacities should hopefully 
make them available in a close future. 



9 IKuchner & S tark (2010) results regarding the respective balance 
between PR drag and collisions are more complex than this simplified 
statement, and depends on the value of the system's average r as com- 
pared to a critical "crossover" optical depth r r 
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5. Conclusions 

Using the DyCoSS code, we numerically investigate the effect of 
collisions on the formation and survival of planet-induced struc- 
tures in debris discs. We confirm the qualitative result found by 
some previous studies using less sophisticated models, i.e., that 
collisions have a global tendency to even out sharp dynamical 
structures. We are here, however, for the first tima^ able to take 
this general result a step further by quantifying this effect and its 
consequences on observability. 

For a planet embedded in an extended disc, we find that the 
planet's signature in the disc could probably remain visible if 
the system is seen head-on, provided that the planet is either on a 
nearly circular orbit and has a mass > 2 x 10~ 4 M*, or that it has a 
mass > 2x 10~ 3 M* if it has a e p ~ 0.2 orbit. If the system is seen 
edge-on, the situation is much less favourable: for bright and 
collisionally active discs, collisions are able to erase most small- 
scale signatures induced by the presence of a planet, although a 
small hollow in the luminosity profile might be detectable un- 
der very favourable conditions. Large-scale structures, in partic- 
ular a steepening of the radial luminosity profiles interior to the 
planet's orbit, should in principle be detectable, but they might 
be very difficult to unambiguously associate to the presence of a 
planet, as other causes could mimic this effect. 

For the other configuration we considered, a planet exterior 
to an inner ring, we find that the planet is never able to prevent 
collisions and radiation pressure from populating the regions be- 
yond the main ring with small fragments, even when these re- 
gions are dynamically unstable. In fact, the luminosity profile 
beyond the main ring is in most cases relatively close to what 
it should be in a planet-less case, i.e., decreasing as r -3,5 . Only 
for tenuous discs with optical depths in the < 2 x 10~ 4 range 
do we obtain significantly steeper profiles. There is, however, a 
narrow transition region just outside the main ring where steep 
luminosity profiles can be obtained, but here again the presence 
of a planet does only moderately affect their slope. The planet 
can, however, affect the morphology of the region correspond- 
ing to the main parent body ring. In this region, the signatures 
of a planet on a circular orbit are azimutal inhomogeneities that 
precess with the planet. These features remain visible despite the 
blurring due to the coupling of collisions and radiation pressure. 
Detecting and measuring the precessing rate of such features in 
observed discs could be an unambiguous way to indirectly in- 
fer the presence of an outer planet, and even constrain its lo- 
cation. If the planet is eccentric, then its signature is to render 
the main ring elliptic and to induce the well known pericentre 
glow effect (periastron side brighter than the apoastron side), 
which is still clearly visible even in highly collisionally active 
systems. However, such features, if detected in real discs, can- 
not be unambiguously attributed to a putative outer planet, as 
other causes, such as an inner planet close to the "natural" outer 
limit of the planetesimal disc, could induce similar signatures. 

For both considered set-ups, lowering the collisional activity 
in the disc to values of to of a few 10~ 4 results in more pro- 
nounced spatial structures. In other words, for a given planet 
mass and orbit, there is a trade-off between the brightness of the 
disc and the sharpness of planet-induced structures: the brighter 
the disc, the higher the collisional activity and thus the lower the 
sharpness of the spatial features. 

DyCoSS, together with the CGA o flStark & Kuchnerl(l2009l) . 
is a first step towards the next generation of debris discs mod- 

1 With the notable exce ption of the collisional grooming algorithm 
of Stark & Kuchner ( 2009), which is, however, intended for a different 
set-up and has been used for a different range of planet masses 



els, where collisions and dynamics will be self-consistently in- 
tegrated into the model's structure. 
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